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Abstract — It has been shown through simulation results 
[14]-[16] that the interference adaptive dynamic channel alloca- 
tion (IA-DCA) scheme is a promising resource allocation strategy 
in time/frequency-division multiple-access (TDMA/FDMA) com- 
munication systems. The major obstacle in analyzing IA-DCA is 
the computation of cochannel interference without the constraint 
of conventional channel reuse factors. To overcome this difficulty, 
one needs a computationally efficient representation which can 
approximate the interference distribution accurately. For this 
purpose, a concept called channel reuse zone (CRZ) is introduced. 
Based on this new concept, both downlink and uplink cochannnel 
interference are computed with two different propagation models, 
namely, a simplified deterministic model and a shadowing model. 
The results are then used to calculate the outage probability of 
the idealized, interference adaptive maximum packing (IAMP) 
scheme. Finally, as a significant contribution, an asymptotic 
performance bound for the two-way IA-DCA strategy is derived. 

Index Terms — Channel assignment, cochannel interference, 
resource allocation, wireless networking. 

I. Introduction 

ONE OF the most difficult problems in the wireless 
communications area is the scarcity of radio resources. 
In order to meet the challenge of the rapidly increasing de- 
mand for the system capacity researchers are exploring new 
resource-saving techniques in all layers. Channel allocation is a 
critical task performed by the network layer in bandwidth-lim- 
ited networks, i.e., time/frequency-division multiple-access 
(TDMA/FDMA) systems. 

In most existing cellular networks, a fixed number of chan- 
nels is assigned to each cell. The same channel cannot be reused 
within a reuse distance in order to guarantee the carrier-to-in- 
terference ratio (CIR) of a connection between a base station 
(BS)/mobile station (MS) pair. This is called the fixed channel 
assignment (FCA). Various dynamic channel allocation (DCA) 
strategies, which aim to enhance the efficiency of bandwidth 
usage, have been devised since the idea was first suggested in 
early 1970's [1], [2]. In the basic DCA, any channel can be 
assigned to any cell as long as the CIR constraint is satisfied. 
Furthermore, other network dynamics can be exploited. 
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The existing DCA strategies are classified into three cate- 
gories based on the type of network dynamics they exploit and 
the representation of the CIR constraint they employ. Traffic 
adaptive DCA (TA-DCA) strategies [3]-[7] take advantage of 
the unevenness of traffic volume among cells and use concepts 
from the first-generation FDMA networks, such as "channel 
reuse factor" and "compatibility matrix," to represent the CIR 
constraint. Another type of DCA which exploits the mobility 
of users is location adaptive DCA (LA-DCA) strategy [8], 
[9]. In the typical implementation of LA-DCA, i.e., reuse 
partitioning, a cell is split into a number of concentric subcells 
and the CIR constraint is approximated by different channel 
reuse factors for different subcell groups. The third type of 
DCA which assigns channels based on real-time interference 
measurement is called interference adaptive DCA (IA-DCA) 
strategies [10]— [13]. There is no formal representation of the 
CIR constraint in IA-DCA schemes. More comprehensive 
overviews on DCA can be found in [17]— [19]. 

Since the received CIR is the ultimate constraint in the 
resource allocation problem, it is reasonable to predict that 
IA-DCA methods will achieve more capacity gain than other 
channel allocation techniques. In fact, many simulation studies 
have confirmed this prediction [13]— [16]. However, theoretical 
analysis of IA-DCA has not been thoroughly addressed in 
the literature. The idea of approaching the upper bound of 
the system capacity by dividing a cell into infinite number 
of concentric zones has appeared in the literature [20]-[22]. 
However, since the method of calculating the perimeter of an 
arbitrary zone for a channel reuse factor other than the classical 
values (1, 3, 4, 7, . . .) was unknown, the authors in these refer- 
ences could only assume hexagonal or circular zone shape and 
could not alleviate the limitation of the conventional channel 
reuse factor. This means that the bounds they obtained were 
for LA-DCA rather than IA-DCA schemes. The performance 
bound of IA-DCA in one-dimensional (1-D) highway networks 
was investigated in [23] (without power control) and [24] (with 
power control). However, it is hard to extend this decent work 
to the two-dimensional (2-D) planar case. 

The main obstacle for the performance analysis of IA-DCA 
in planar networks is the representation of cochannel inter- 
ference (CCI) constraint without approximation. The existing 
literature on CCI computation addresses the code-division 
multiple-access (CDMA) systems [25]-[27] and FCA-based 
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TDMA/FDMA: system [20], [27]. In this paper, we present a 
method for computing the concentric zone boundaries corres- 
ponding to any nontrivial value of channel reuse factor.' The 
result is: not only are we able to get rid of the confine of the 
conventional concept of the channel reuse factor, but also we 
can extend our investigation from the downlink-only DCA, 
deterministic propagation model (without shadowing) to the 
case of both the link and shadow fading models. 

The contribution of this paper can be summarized as fol- 
lows. First, a novel concept called the channel reuse zone (CRZ) 
and its extension to the shadowing model, the extended CRZ 
(ECRZ), are introduced. Furthermore, a CRZ/ECRZ structure 
called the concentric CRZ structure is described. This struc- 
ture represents the most compact channel reuse pattern under 
the CIR constraint. Second, an idealized IA-DCA scheme called 
interference adaptive maximum packing (IAMP) algorithm is 
defined. The motivation of the study on IAMP is to find the per- 
formance bound of the IA-DCA strategies. The CCI and outage 
probability are calculated for IAMP in both downlink and uplink 
channels. Finally, a performance measure called the asymtotic 
probability of assignment failure, introduced in [20], is used as a 
benchmark and its lower bound (the upper performance bound) 
for IA-DCA strategies is derived by applying results obtained 
in the first two parts. All of the computations and derivations 
consider both the deterministic and shadowing models. 

The paper is organized as follows. Notation, assumptions and 
propagation models are presented in Section II. Section III fo- 
cuses on the deterministic model. First, the CRZ concept and the 
IAMP algorithm are given, and then the computation of CCI is 
presented. The results for the deterministic case are extended 
for the shadowing model in Section IV. Section V deals with 
the derivation of the lower bound of the asymtotic probability 
of assignment failure for IT-DCA. strategies. Numerical results 
and interpretations are presented in Section VI. 

II. Notation, Models, and Assumptions 

The typical cellular network architecture is considered in this 
paper. A cell is a geographic area corresponding to the radio- 
service area of a BS. 2 An MS communicates with an assigned 
BS through a channel with an assigned transmit power level. 
A channel might be a frequency band in FDMA systems or a 
time slot in TDMA systems. The effect of power control is not 
considered in this paper. 

The CIR constraint rules that a channel p can be assigned to 
a BS/MS link only if (C/I) p < 7, where (C/I) p is the CIR 
of channel p and 7 is the desired CIR threshold. A downlink 
(uplink)-only DCA algorithm assigns a channel pair if the CIR 
on the downlink (uplink) channel is above the 7. This algo- 
rithm assumes that the CIR's on both channels are identical. A 
two-way DCA scheme assigns a pair of channels only if both 
downlink and uplink channels meet the CIR constraint. 


'A "nontrivial value" means a real number in the range of [1, h K ] , where h K 
is determined by the CIR threshold and the propagation model (see definition 
in Section III). 

2 A11 the theoretical results in this paper are independent of the cell shape. For 
the sake of simplicity, the numerical results in deterministic case (in Section VI) 
are acquired by assuming hexagonal cells. 


A. Notation 

The following notation is used in the paper. 
R Radius of a cell. 

A c Area of a cell. 

D Distance between two neighboring BS's. 

D = VdR. 
rji Distance between BS J and MS i. 

P d (P u ) Received downlink (uplink) local mean signal 

power. 
Lji Shadowing effect of the received signal. 

a Path-loss exponent. 

Uj Number of ongoing calls in cell j. 

F Average rate of blocked and dropped calls. 

M Total number of channels in the system. 

V Number of cells in the system. 

P out Outage probability. 

B. Propagation Model 

The downlink signal power received at MS i from BS J is 
expressed as 


p ji - PoLjir, H 


—a 


(1) 


where Pq is the received signal power at unit distance. Similarly, 
the uplink signal power received at BS I from MS j is given as 


Pfr^PoLjjrTf. 


(2) 


In our study, we assume that: 1) the system is interference lim- 
ited, such that back noise power is negligible; 2) the shadow 
fading Lji's are independent lognormal random variables (RV) 
with a mean of dB and a standard deviation of a dB 3 ; 3) omni- 
directional antennas are used; and 4) interchannel interference 
is negligible, so only cochannel interference is considered. 

C. Traffic Model 

We assume that the traffic in the system has uniform distribu- 
tion, hence Uj is a random variable with statistics independent 
of cell index j. Furthermore, offered traffic in each cell is as- 
sumed to be an independent stochastic process with an average 
arrival rate A (calls/cell/ s). The mean duration of each call is 
1/p, (s). The average offered traffic in a cell is 


£ = A//u.(calls/cell). 


(3) 


We can define the relative traffic load as p = £,/M [20]. The 
probability of assignment failure is defined as v — E[F]/£V. 
Therefore, at any snapshot in time, Uj's are independent RV 
with a mean of A(l — v)/ p,. 

A keen reader may notice that the mobility of mobiles is not 
explicitly considered in this model. In most cases, the mobility 
issue makes a huge difference in the design and evaluation of a 
practical DCA algorithm. Mobility increases the signaling load 
in the system. It may also lead the DCA designer to give pri- 
ority to handoff calls at the cost of reducing the system ca- 
pacity. Nevertheless, mobility is not explicitly modeled in this 
study due to the following reasons. First, the impact of mobility 

'The i.i.d. assumption is made for the sake of simplicity. The result can be 
easily modified for the correlated case by applying formulas in [28] . 


XU et al.: ON DEPLOYMENT REQUIREMENTS OF OVERLAID LTE TRANSCEIV- 


783 


on system capacity is represented by increased rate of channel 
seizure attempts (A) accompanied by shorter channel holding 
times (1//U,). The resulting snapshot traffic load £ in (3) is similar 
to £ without mobility. Second, in order to find the performance 
bound for IA-DCA strategies, one needs an allocation algorithm 
which is more powerful than any other IA-DCA schemes. For 
this purpose, we design the idealized IAMP method, which is 
assumed to have unlimited processing power to perform all of 
the necessary channel allocations/reallocations instantaneously. 
Therefore, the probability of handoff failure due to processing 
delay is zero. Third, the probability of assignment failure v, 
the benchmark performance measure in our study, includes both 
probability of new call blocking and forced termination and does 
not distinguish between them. 

To calculate the performance measures for IA-DCA schemes, 
one needs to know the maximum number of users the system 
can accommodate. For this purpose: 1) a representation of the 
CIR constraint has to be defined in an IA-DCA environment 
and 2) the outage probability at every point in the cell has to be 
calculated for the given DCA algorithm and the given represen- 
tation of the CIR constraint. The calculation of P out depends 
on the CCI computation. The next two sections are organized as 
follows: 1) the CIR constraint; 2) the representation of the con- 
straint and definition of the allocation algorithm (IAMP); and 
3) the computation of interference and P out for both uplink and 
downlink. 

iii. cochannel interference analysis for 
Deterministic Case 

In the deterministic case, shadow fading in (1) and (2) is set to 
be dB. The downlink and uplink CIR constraints for the cell 
of interest (denoted by cell 0) are thereby defined as 


{C/If 


'oo 


>7 


(4) 


and 


Jei d 


'JO 


(c/iy 


'00 


respectively, where 




>7 


(5) 


cells sharing the downlink channel with cell 
and l<i = { 0} 

cells sharing the uplink channel with cell 

J" = { 0} . 

A. Representation of the CIR Constraint — Channel 
Reuse Zones 

Imagine the scenario that a mobile on channel moves from 
BS to the border of the cell in an arbitrary direction. At the 
beginning, the MS is very close to BS so that (9) holds 
even if all other cells are transmitting on channel 

p. When MS moves 
away from BS 0, its CIR decreases. At a certain point, one of the 
neighboring cells of cell has to be forbidden from transmitting 


on channel p, in order to satisfy (4). On its way to the cell border, 
MS causes more and more cells to be prohibited from using 
channel p. Since the BS's are discretely located, CRZ's are nat- 
urally formed. 

Definition 1: A downlink CRZ (CRZ) is a region in a cell 
where (4) holds for the same interference distribution despite 
of the mobile 's location. An uplink CRZ is a region in a cell 
where (5) holds for the same interference distribution despite 
of the mobile's location. A two-way CRZ is a region in which 
both (4) and (5) are simultaneously satisfied. 

It is clear that the shape and area of a CRZ are determined by 
a, 7, and the distribution of cochannel interferers. The number 
of possible zone partitionings is virtually infinite. A special par- 
titioning of CRZ's is defined below. It turns out, by CCI com- 
putation (see Section III-B) that this specific CRZ partitioning 
consists of a set of concentric zones. Therefore, we call it as 
a concentric CRZ structure. The word zone in the next para- 
graph stands for either a downlink, uplink, or two-way CRZ. 

In the concentric CRZ structure, zone-1 is defined as the 
largest CRZ in which a channel can be reused in all the cells 
without violating the CIR constraint. Then, the cells in the 
system are divided into V/2 groups. Each group has two 
neighboring cells. Zone-2 is defined as the largest CRZ in 
which a channel can be reused in one of the two cells in a cell 
group. The word "largest" means that the cochannel cells from 
different groups should be distributed as sparsely as possible. 
By dividing the system into V/3 identical groups, with three 
neighboring cells in each group, we can define two CRZ'S: 
Zone-(3/2) is the largest CRZ in which a channel can be reused 
in two out of three cells in a cell group. Zone-3 is the largest 
CRZ in which a channel can be reused in one out of three cells in 
a cell group. In general, Zone-(n, 1) (n, I integers, n > I > 0) 
is obtained by splitting the system into V/n identical cell 
groups, each with n cells, and assigning a channel to I cells 
in each group in such a way that the resulting CRZ covers the 
largest possible region (or the cochannel cells are distributed 
most sparsely). Let K be the number of CRZ's in each cell, 
the Kth zone is denoted as zone-hk, k — 1, 2, • • • , K, where 
hi - 1, h k > h k -i- 

The physical meaning of hk is that a channel can be reused in 
one out of hk cells within zone-hk ■ Therefore, it can be viewed 
as an extension of the "channel reuse factor." However, it should 
be emphasized that hk might take improper fraction values. This 
implies that K may go to infinity, although the value of hk is 
upper bounded by a and 7. Note that for group indexes n > 2, 
the layout of the cell group is not unique. For example, four cells 
may be arranged in a row or stacked in two rows. Even the shape 
of a cell group which gives the largest CRZ area may change for 
the same k. However, for a given hk, the area of the largest CRZ 
is identical in every cell since the interference distributions are 
the same for every cell group. Since only the area of CRZ'S, 
which corresponds to the offered number of calls, is useful in 
our analysis, the definition of zone-hk is reasonable. 

Figs. 1 and 2 illustrate the relationship between zones in the 
concentric CRZ structure and the corresponding channel reuse 
patterns. Numbers in cells represent the distance between the 
cell of interest and another BS (unit: D). Fig. 1 shows the most 
sparse cochannel user distribution for a channel used in every 
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Fig. 1. Channel reuse pattern of concentric CRZ zone-5. 



Fig. 2. Channel reuse pattern of concentric CRZ zone-13/2. 

one out of five cells. This reuse pattern is utilized to compute 
the area of CRZ zone-5. Fig. 2 shows the most sparse channel 
reuse pattern corresponding to zone-13/2. 

Definition 2: A zone-hk user is defined as a mobile located 
within zone-hk and out of zone-hk— i- A zone-hk channel is a 
channel assigned to a zone-hk user, which might be either an 
uplink or a downlink channel. 

The theorem below is the basis for the study on the perfor- 
mance bound of IA-DCA. 

Theroem 1: Assume that: 1) traffic is uniformly distributed 
in each cell and 2) the allocation policy is fair for all cells, i.e., 
the blocking probability in each cell is the same for equal traffic 
loads. As the number of CRZ' s K goes to infinity, the concentric 
CRZ structure represents the CIR constraint of the most compact 
channel reuse pattern in TDMA/FDMA systems. 

Proof: The most compact channel usage in TDMA/ 
FDMA systems is the one where a channel is reused in all V 
cells (reuse pattern 1). To use channels in an optimum way, 
reuse pattern 1 should contain as many users as possible. 
In other words, reuse pattern 1 should cover the averagely 
largest possible region in a cell (by assumption a), and the 
area of this largest region should be identical in every cell (by 


assumption b). This argument results in the definition of the 
concentric CRZ zone- 1 . For the users outside zone- 1 , the next 
most compact reuse pattern is when one channel is reused in 
(n — l)V/n cells, where n is an arbitrarily large positive integer. 
By the same argument as before, we obtain the definition of 
zone-n/(n — 1). By continuing the same procedure on users 
outside zone-n/(n — 1), one obtains the whole concentric CRZ 
structure in a straightforward way. 

It is apparent that the concentric CRZ structure with finite K 
is still an approximation to the CIR constraint. However, this 
structure can approximate the original constraint with an un- 
limited accuracy by increasing K. The real CIR constraint of 
the most compact channel reuse pattern will be reached when 
K — > oo. 

B. The Interference Adaptive Maximum Packing Method 

In order to study the performance bound for IA-DCA 
schemes, every channel is expected to be packed as densely 
as the CIR constraint allows. Due to the limited processing 
capability of a real network, no practical DCA scheme is able 
to take full advantage of the channel reuse pattern represented 
by the concentric CRZ structure. Therefore, the idea of the 
well-known maximum packing (MP) approach (see [5]) is 
employed for the purpose of theoretical study. The basic idea 
of MP is of assuming that there exists a central controller who 
knows all the system-wide information and has an unlimited 
computational power. The controller will do all of the necessary 
channel reallocations in the whole system in order to accom- 
modate a new call. Since MP is originally used in TA-DCA, 
the method proposed here is called the IAMP method to avoid 
confusion. The following algorithm may be viewed as either 
downlink/uplink-only or two-way IA-DCA, depending on what 
type of CRZ zones is used. 
Step 1) Set k = 1. 

Step 2) Pick one zone-hk user from every group of hk 
cells, do all the necessary channel reallocation, and 
pack them into a single channel (or a channel pair 
in two-way case). The cells occupying the same 
channel should be in the same location within their 
own groups (see Figs. 1 and 2). In the two-way 
case, uplink and downlink channels are assigned 
such that cells using the same uplink channels are 
assigned the same downlink channels and vice 
versa. Continue this process until all zone-hk users 
are packed or all available channels are used up. 
Step 3) Increase k by one; go back to step 2 if k < K, 

otherwise, stop the algorithm. 
The idea of ordering the channel packing is to increase ef- 
ficiency. Since a zone-hk channel can only be reused in V/hk 
cells, the maximum number of users a channel can support de- 
creases as k increases. Therefore, it is more efficient if channels 
are assigned to zone-hk users first, then zone-/i2 users, zone-h^ 
users, and so on. 

The significance of binding uplink and downlink channels 
together can be conceived through the following inequality: 


i^ t =l-(l-i? ut )(l-i* 

>max[Ct,^t]=^ 


out/ 

bd 
out 
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Fig. 3. Computation of downlink cochannel interference. 

where F" ut , P^ ut , i^„ t , and.P "^ t are outage probabilities of the 
uplink channel, downlink channel, two-way channels assigned 
with binding, and two-way channels assigned without binding, 
respectively. 

Theorem 2: Under the assumptions in Theorem 1, the IAMP 
method applied to the concentric CRZ structure becomes op- 
timum (in terms of channel use efficiency) in DCA strategies as 
K goes to infinity. 

Proof: From Theorem 1, when K goes to infinity, the 
concentric CRZ structure represents the most compact channel 
reuse pattern. In other words, it gives the largest average set of 
available channels under the CIR constraint. The assumption 
of unlimited processing power of IAMP implies that if there is 
an available channel in the system, the IAMP method is able to 
assign it within an arbitrarily short time period. As the result, 
the combination of these two factors gives the optimal channel 
utilization. 

C. Computation of Cochannel Interference: 
Deterministic Model 

I) Computation of CIR for General Case: The computation 
of downlink cochannel interference is explained through Fig. 3. 
To simplify the notation, let us define r — roo and 9 and <f>j 
as shown in Fig. 3. The downlink local mean interference is 
expressed as 


Pjo - For jo 


-(a/2) 


= Po[w]D 2 + r 2 -2w j Drcos(4> j -9)] v * / "' (6) 

where Wj = 1, a/3, \/7, 3, • • •, depending on the distance be- 
tween BS J and BS 0. The downlink CIR at MS is written as 


(CIR) d 


r 


Y^ [w 2 p 2 + r 2 - 2w j Drcos(<f> j - 9)] 

)El d 


-{a/2) 


(7) 

The computation of the uplink interference is illustrated in 
Fig. 4. Now, the difficulty is that every vjo is an independent 



Fig. 4. Computation of uplink cochannel interference. 

RV. Therefore, an equation like (6) which relates PJ with the 
position of MS does not exist. The main difference between 
the uplink interference computation in CDMA (see, e.g., [25] 
and [26]) and this study is that in CDMA the interference from 
another cell is due to an accumulated effect of average £ users, 
whereas in TDMA/FDMA there is up to one interferer in each 
cell. Therefore, instead of computing the expected total interfer- 
ence from a cell, we compute the expected interference power 
received at BS from MS j as 


P* =E{Portf} 


= PoE 


| (w 2 D 2 + r 2 - 2wjDtj cos < 


-(«/2)j 


(8) 


where the expectation is over a specific region in which MS 
j may appear. In cases of FCA and traffic adaptive DCA, this 
region covers a whole cell. The uplink CIR at BS is expressed 

as 

,»,— a 

(CIR)" 


£*{( 


w 2 } D 2 + rj 


— 2wjDrj cos ( 


-(a/2)' 


jei u 


(9) 
2) Computation of CIR and P ou t for IAMP Method: 

a) Downlink CRZ: Under the IAMP method, the border 
of a downlink CRZ zone-hk can be found by solving the radius 
RkdiP) f° r various (?'s from the following: 


R 


kd 




Jeii 


[w]D 2 ■ 


R 2 d -2w j Drcos^ j -9)] 


-(a/2) 


(10) 

This equation is obtained by: 1) defining the most sparse 
cochannel user distribution corresponding to hk', 2) inserting 
(7) into (4); and 3) setting the inequality in (4) to equality 

d emphasizes that 


and setting r — Rkd- The subscript k of X h 


it changes with different hk- The concentric CRZ structure 
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uplink zones 


downlink zones 


Fig. 5. Concentric CRZ structure with integer hk's. 

is formed due to the fact that interferers move farther as hk 
increases. 

b) Uplink CRZ: In the case of IAMP, all of the uplink 
cochannel interference in a zone-hk channel is from zone-hk 
users in other cells. As a result, the expectation in (8) becomes 


P]=P 


2?r 


Rk 


(w 2 D 2 + r 2 - 2wiDn cos < 


-(a/2) 


'0 JR k - X 

■Pk(r 3 ; fa)dr 3 -dfa 

p r 2pi r R k 


riRl-RU) Jo Jr^ 


■ (wjD z + rj - 2wDvj cos fa) 


-(a/2) 


rj dr k d(f>j 


(11) 


where pk(r } -, fa) is the probability density function (pdf) of a 
zone-hk user being at (r } -, fa) (uniform by assumption). Then, 
the radius of uplink zone-hk, Rku, can be solved from 


^{Rlu-Rlu- 


E 


Rhu 


jeI uJRku- 

= 7- 


(w 2 D 2 + rj 


2wjDrj cos 4>j) " Tj drj dcftj 


(12) 


Equation (12) is obtained by following a similar procedure as in 
the downlink case. Note that (12) implies that Rk u is indepen- 
dent of 9 and thus uplink CRZ's are circles. 

Fig. 5 shows the concentric CRZ structure of the IAMP 
method in the deterministic case, where hk's are integers and 
a = 4.0 and 7 = 18 dB . By definition, a two-way zone-hk 
is the overlap region between downlink zone-hk an d uplink 
zone-hk- 

The computation of P out is rather straightforward in the CRZ 
structure. For a specific interference distribution corresponding 
to hk, the outage probability P out equals zero within zone-hk 
and equals one outside of it. 


IV. Cochannel Interference Analysis for 
Shadowing Environment 

In the shadowing environment, the downlink and uplink CIR 
constraints are expressed as 

L r~'- 


{c/iy = 


> 


E^ 

Jei d 


1 


jor 


JO 


and 


{c/iy 


L r' 


E 


Ljorj? 


>7 


(13) 


(14) 


respectively. In the deterministic case, the location of a mobile 
is the sole factor that determines its CIR for the given 

a and 
7 parameters. This factor is mapped into CRZ's. In this way, 
the CIR of a zone-hk user is guaranteed by assigning a zone-hk 
channel to it. The only reason left to block a call is the lack of 
available channels. However, in the shadowing model P out is 
no longer only zero or one as in the deterministic model. It is a 
function of hk, r, and 9. As a result, the relationship between 
the geometric borders (and the area) of CRZ's and the number 
of zone-hk users is completely destroyed. As we will see later, 
this relationship plays a pivotal role in the study of performance 
bound. In order to recover this relationship, the concepts of CRZ 
and cell have to be revamped as follows. The new definitions 
call for the terms "the extended cell" and "the extended CRZ" 
for consistency of terminology, although there no longer exist 
"zones" in the geographic sense. 

A. Extended Cell and Extended CRZ 

Regardless of whether there is shadow fading or not, an ac- 
tive MS always has a higher CIR with one of BS's at a given 
moment. Since the essence of concepts CRZ and cell is to rep- 
resent the CIR constraint rather than the geometric distance, the 
definitions should be modified according to (13) and (14). 

Definition 3: A downlink (uplink) extended cell (EC) ofBS 
J is a group of active mobiles, with a mean value £, which has a 
higher downlink (uplink) CIR on a channel with BS J than with 
any other BS. 

Definition 4: A downlink (uplink) ECRZ ofhk, denoted as 
zone-hk, is the largest group of mobiles in a downlink (uplink) 
EC whose CIR > 7 on a channel assigned to one out ofhk cells. 
A two-way ECRZ ofhk is defined as a group of active mobiles 
which belong to both the downlink and uplink ECRZ zone- hk- 

The word "largest" means that if there exist more than one 
interference distribution patterns for a specific hk, then zone-hk 
corresponds to the one with the lowest P out curve. 

Obviously, it is extremely difficult to identify EC's 
and ECRZ's one by one, even for a snapshot. However, 
it is not that difficult to find the "area" of an EC or an 
ECRZ. The "area" of an EC/ECRZ is defined as the av- 
erage number of users in the EC/ECRZ at a given time 
instance. For the purpose of simplicity, let us assume that 
the solely function of the distance r. Thus, (the downlink) 
P out (r) = F r {the CIR of an MS at distancer < 7}. Now the 
question is: if we put m users at the distance r from a given 
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BS, what is the average number of successful users who is 
not blocked? It is straightforward to show that the average 
number of successful users at distance r — X)n=o n \)- ~ -Pout]- 
Therefore, the "area" of ECRZ zone-Zi^ equals 

/ / r[l-P out (h k ,r,e)]drde (15) 

Jo Jo 

where P ou t(hk, r, 9) depends upon the specific DCA algo- 
rithm. Rfr is the upper bound of r, such that P ou t(hk, Rk) ~ 1- 
The "area" of an EC is the summation of "areas" of ECRZ's 
from one to K. 

B. Computation of Cochannel Interference: Shadowing Case 

1) CIR Computation for General Case: 

a) Downlink CIR: From (6) and (13), the downlink CIR 
at user in shadowing case is express as 

(CIR) d 

L r~ a 


where Vd is the number of cells in I d . After fj,^ and [i 2 d are 
obtained, m z <i and a 2 d can be found as 

m zd = 2 In ma - \ In fi 2 d (22) 

o- 2 zd - m M2d - 2 In /j, ld . (23) 

Since both Y and Zd are Gaussian, Y — Zd is also Gaussian 
with 

N f -m zd ~ alnr, y a 2 + a\ d \ . 

Then, the downlink outage probability of an MS at /y @\ is 
expressed as 


P* ut (r,e)=P r {Y-Z d <ln 1 } 


-l-O A I1 T + alnr + m «d 


(24) 


In (24), can be represented aya^funcS-efh of / . As a 


result, 


and 


are functions of 


iftheinter- 


Y^ Ljo [w]D 2 + r 2 - 2w j Drcos(<f> j - 9)] 

)El d 


-{a/2) 


■pfrf 


(16) 
The major problem in computing P^ ut with a shadowing 


ferencedistribution aradcihe parameters / a\ are known. 

b) flplinttlR: From (8) and (M), the Vplink CIR at BS 
is 


model is dealing with Ljq's. Let 


= L r- a e Zj = L 


CIR 

( )'" 


J0' r J0 


vj eT 


d 


(17) 


Lor~ a 


Since L } 's are assumed to be lognormal RV's, Y and Zj are 
both Gaussian with N(m y , a 2 ) and N(mj, a 2 ), respectively, 
where m y — —alnr and nij — — alnrjo. Moreover, let e Zd 
be the summation of lognormal RV's as 


*/2)\ 


Y L }0 E{{w]D 2 + r]-2w } Dr } - cos <f> } ) {a/ j 

(25) 


jei u 


,-,z d _ y^ „z,, 

}<El d 


(18) 


Combining (16)— (18), we get 

^out = Pr {e Y /e Zd < 7} = P r {Y-Z d < ln 7 } . ( 19 > 
Note that e z d is a weighted summation of a group of log- 
normal RV's. There exist three major approaches to 
approximating such a summation [28], namely, the Wilkin- 
son's approach, Schwartz- Yeh approach, and cumulant 
matching approach. We use Wilkinson's method in our 
investigation due to the conclusion of the comparative study 
in [28] and [29]. In Wilkinson's method, is assumed to 

be another lognorma^RV. In other words, is^jaussian 
with m£an and variance , which are obtained by 


where the expectation is taken over the region in which MS j 
may appear. The region depends upon the specific DCA 
algorithm. 

Assuming that 


rm 


E \ (wW z + ri - 2w i Dr i cos 


»,)-«■'«} 


!-(!/«) 

(26) 

is known for a given DCA algorithm, pu can be found in the 
same way as in the downlink case with obvious changes on 
indexes 


Pout(r, fy ) = 1 - Q 


where 


In 7 + a In r + m zu 


■ zd 


equating first two moments of 
(18) as 


v d 


IJ'ld 


-E{e Zd ) = E \ 5>*')=^ /2 »Jo a 
^Jei d 


/ , 'JO 
J = \ 


(20) 


and 


m zu — 2 In /u," - 0.5 In /j% 


and 


IHd 


--E[e 2Zd ] 



Mi" = e 




(27) 

(28) 
(29) 

(30) 
(31) 


j=i 


^ e 2—i 2—/ ' Jo ' so 

J=l S=J+1 


r^rcn" (21) 


„2<H 


K = e~ E^ + 2e^X: £ rj«r 7 « 

j=l j=l 8=3+1 


788 


IEEE TRANSACTIONS ON VEHICULAR TECHNOLOGY, VOL. 58, NO. 9, MAY 2012 


After both P^ ut and P™ ut are calculated, the two-way outage 


probability is computed after the definition 

Pout(r, 9, f j0 ) = F d ut (r, 9) U i^ t (r, f j0 ). 


(32) 


2) CIR Computation for IAMP Method: With the IAMP 
method, the interference on every downlink/uplink zone-Zi^ 
channel is from other downlink/uplink zone-hk channels, and 
downlink and uplink channels are bound together. Thus, for a 
zone-hk user m ce ll 0, if. = 1% = Ik with interferer number 

A 

V<i — V u — Vk — V/hk — 1. To emphasize the fact that in 
IAMP method the outage probability depends on hk, we write 
F d ut as P« ut (h k , r, 9) and F« ut as P- ut (h k , r, f j0 ) below. 

a) Downlink CIR: With the new notation, P^ ut of a 
zone-hk user at (r, 9) can be found from (24) with ij,^ and 
\i2d changed to 


IJ'id - e 


v k 

2/2 £ 


JO 


and 


J=i 


2ef 


v k v k 

'E E 

J=l S=J+l 


JO 'so 


(33) 


(34) 


respectively. 

b) Uplink CIR: In the uplink case, the expected distance 
between BS and a zone-^ user in ^ . {$ calculated from 

(26) by taking the expectation over ECRZ zone-hk m ce U j 
(assume uniformly distributed traffic) as shown in (35), given 
at the bottom of the page. Note that Rk is the upper bound for 
r 3 such that P^ ut (h k , r 3 ) w 1. 

Equations (26)-(31) and (35) form a family of equations with 
unknown function P™ ut . Rk in (35) is unknown a priori. One 
way to solve this equation is to search for the solution recur- 
sively. A good guess for the initial P™ ut is the corresponding 
Pp Ut found in (24). The numerical experiment shows that P" ut 
can be reached within five-six iterations if P^ ut is chosen as the 
initial function. 4 

After both £ ut and P™ ut for IAMP are obtained, the two-way 
outage probability is computed as the following due to the bound 
downlink/uplink channels in IAMP 

Pout(h k , r, 9, Fyo) 

= max {P^ ut (hk, r, 9), P^ ut (h k , r 3 , f j0 )} . (36) 

4 The stop condition is the maximum absolute error between two consecutive 
iterations is less than If) -4 . 



Fig. 6. Probability of outage under shadowing (symmetric interference). 
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Fig. 7. Effect of 


r/R for various hk with the same a, 7 , 
and a — 8 dB. Fig. 6 shows both P^ ut and P™ ut for hk's whose 
interference distribution is symmetric and therefore 9 has no in- 
fluence on Pout- Note that P™ ut is higher than % ut in most cases. 
To appreciate the effect of shadowing, notice that the determin- 
istic case, with a = 0, zone-1, zone-3, and zone-7 are circles 
with normalized radii of 0.357, 0.622, and 0.951, respectively. 


rjo(h k ) 



-(a/2) . 


-(l/«) 


(w 2 jD 2 + r| - 2wjDtj cos <j>j) Pk(r 3 , </>y) dr 3 d(f>j 


"(«/2) r™ 


w]D 2 + r] - 2w 3 Dr 3 cos cf> 3 ) w "> [i? ut (h k -i , 77) - i? ut (h k , r 3 )] r 3 dr 3 dcf> 3 
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(35) 
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From Fig. 3, it is obvious that because of shadow fading of the 
signal, a mobile's CIR is no longer guaranteed. For example, a 
zone-1 user at a distance of 0.3.R to its BS has an outage prob- 
ability around 0.5. On the other hand, because of the shadow 
fading of the interference strength, some users at distances of 
0.9 ~ 1.0-R might have a chance to be packed into zone-1 chan- 
nels. Such space spreading of users in the same CRZ zone in- 
creases as hk increases. An interesting task here is to check if 
we gain or lose capacity due to the shadow fading. The change 
of expected number of users in each zone will be calculated in 
the next section. 

For CRZ's with hk other than classic channel reuse factors 
(1, 3, 4, 7, etc.), cochannel interference is not symmetrically dis- 
tributed. Thus, mobiles with different angles to the BS have dif- 
ferent Pp Ut . The effect of 9 on the downlink outage probability 
for asymmetric interference distributions is shown in Fig. 7. The 
figure shows that the influence of 9 decreases for either less r, 
because of higher CIR, or larger hk , due to farther interferer 
positions. 


V. An Asymptotic Perrformance Bound for IA-DCA 

As an important application of the CCI computation obtained 
in previous sections, a lower bound of a performance measure 
called "the asymptotic probability of assignment failure" is 
derived for both propagation models. Throughout this section, 
we assume that: $,)M are arbitrarily large, bait is kept 
finite a¥d 2) is large so that the boundary effect is neglected. 
Since we are considering the asymptotic situation, the content 
of this section may be called as asymptotic analysis. 

As we previously mentioned, the key issue in applying the 
concept of CRZ/ECRZ in the performance analysis is to estab- 
lish the relationship between the area of CRZ/ECRZ and the 
number of users it can contain. Definitions in the first section are 
used to establish this relationship. In this section, the word zone 
means either a downlink, an uplink, or a two-way CRZ/ECRZ. 
P ou t means either the downlink, uplink, or two-way probability 
of outage, depending upon what kind of zones are considered. 
Similarly, theorems proven in this section apply to each of three 
IA-DCA methods. 

A. Definitions 

Definition 5: The asymptotic probability of assignment 

failure is defined as v * _ lim ^ u ^^ v 

Definition 6: The area increment of CRZ is defined as 
A k — (area of zone-h k ) - (area of zone- h k -i) 

r 2yr r R k {$) 

r dr d9 (37) 

/o JR k -i(e) 

where Rq = 0. The area change of CRZ is defined as 

A k = A k -i 


a(hk) = lim 


a/u.-»o Ahi 


(38) 


A% = 


Definition 7: The extended area increment of ECRZ is 

defined as 

f 27T r-Rkifi) 

/ [P ut(/ifc-i, r, 9) 

io Jo (39) 

-Pout(h k , r, 9)]rdrd9 
where A 

■Pout {ho j r) — 1 and s stands for "shadowing. " Rk is the 
upper bound ofr, such that P ou t(hk, Rk) ~ 1- The extended 
area change of ECRZ is expressed as 

a 3 (/i fc ) = lim h ~, h ~ 1 . (40) 


Ah k 


B. Lower Bound of for Deterministic Model 
v* 
Theorem 3: The lower bound of the asymptotic probability 

of assignment failure for IA-DCA in a deterministic case can be 

approximately represented as 


fo, 


r/(q - 1) - 


+C(<z), 


1- 


P < IKK) 

l/rj(q) < p 

<iM?-i), 
p > iMi) 


(41) 


where 


v(q)^}~ c t / hk/Ak C(</) = ^X> fc . (42) 

k=l k=q 

Proof: Since Uj is uniformly distributed over a cell 

(Avg. # of CRZ zone-hk users in system) = £VAk/A c . 

(43) 
The maximum number of zone-Zi^ users a channel supports is 
V/hk ■ Hence, the average number of channels needed to support 
zone-hk users is ^hk/Ak/A c . Thus 

K £,h k A k 


total m. of channels needed 


E 


WK). 


Let us define _, the number of failed zone-, users. 

Due to the ordering oFchannel packing in IAMP algorithm, a 


user cannot fail unless all zone- 


users have 


ZOne u»v. wu 1 ,v lll m»u 1 v»» l » lw uv , . 

alreaay^failed. Therefore, when the number of available chan 


d^f 

^S tv(K), < m = o^ And when s v (K) >m> wk - 1) 

E[F\=E[F K ] 

— (zone- hk channel deficit) 

(^ of zom-hk users P er channel) 


^2 h kA-k 


A'-l 


where Ahk — hk — hk — 1. 



Y^ h kAi 


-^Lk-i)- 1 

h K \ P 
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For the case where £rj(K - 1) > M > £rj(K - 2) 


E[F]=E[F K - 1 ] + E[F K ] 


A c 


A'-l 


Y^ h k A-k-M 


V 


k=l 


Hi-l 


■& 


A c 


K-2 


=* v hrE^ = 


M 

J, 


1=ZV 


Ak + M-A 
ac J 

J—L{K-2)-- 
h K -i V p 


h K -i 


<{K-1) 


In general, when £rj(q) > M > £rj(q — 1), q — 2, 3, • • • , K 

E[F] = J2E[FK]=tvt [ ^^(q-l)-^+aq)}- 
Finally, when M < £77(1) 
£[F] = £ £[F A -] 


fc=i 


= ^1- 


1 


A 2 


A A - 


Equation (41) is obtained via the definition of v* . 

Theorem 4: The lower bound of the asymptotic probability 
of assignment failure oflA-DCA in deterministic case is 


— < 


CO/), 


P < IKK) 

1/77(20 < P < 1M1) 

p > lMl) 


(44) 


where 


<(<?) = 

%) = 


1 
^ Jh q 

Ax 


a(x) dx 


xa(x) dx 


(45) 


and hn and h q are defined, respectively, by 


i{x) dx — A c — A\ 


xa(x)dx — — — A\. (46) 


Theorem 5: The lower bound for the asymptotic probability 
of assignment failure of I A- DC A under lognormal shadowing 
can be approximately represented as 


0, 


rj(q - 1) 


+C(3), 


I- 1 , 
P 


P < IKK) 

iKq) <p< l/77(g - 1), 

p > iMi) 


where 


v(q) 


' k=l 


C(<?) 


1 A " 

k=q 


At, 


(47) 


(48) 


■d0 

(49) 


Proof: A tiny region in which P ou t(hk, r, 9) is approxi- 
mately constant for a fixed hk has an area A.A = rArA9. Since 
Uj is uniformly distributed over a cell, the average number of 
active users in A A is £AA/A C . Therefore, average number of 
zone-hk users in the system is expressed as 

- (1 - P out (h k -i, r, 9))] 
£V i' 2lx f Rk 

A / ^t^- 1 ' r ' ^ _ Fo ^ fc ' r ' ^ r dr ( 

c t/ U v U 

The rest of the proof follows a similar way to the one used 
to prove Theorem 3 in deterministic case by making obvious 
changes on Ak and a(hk). 

The final proposition can be obtained by introducing the def- 
inition of a s (hk) and letting K — > 00, as in the deterministic 
case. 

Theorem 6: The lower bound for the asymptotic probability 
of assignment failure of IA-DCA for the shadowing model is 

(0, p< IKK) 

C(q), IKK) >p< 1/77(1) 

1--, p>iKV 


(50) 


where 


<(<?) 


77(g) = 


1 

rhk 

Ac, 

'h q 

1 

A[ 


a s (x) dx 

'■h„ 


} (x) dx 


(51) 


Proof: Using the definition of the area change of CRZ and 
noticing that rj{k) — rj{k — 1) —> as K —> 00, one can obtain 
the above theorem in a straightforward way. 

C. Lower Bound ofv* for Shadowing Model 

The major impact of shadow fading on performance analysis 
is the damage to the simple relationship of (43) in the determin- 
istic case. However, whenever P out is found, the relationship 
is recovered by the introduction of ECRZ as shown in the fol- 
lowing proposition. 


and hn and h„ are determined, respectively, by 


h K r-h q 

a 8 (x) dx — A c - A{ 
1 A 


A c 


Xd \pC ) dX — ■'■M * 

(52) 


VI. Numerical Results and Discussions 

Theorems 3 and 5 give concise expressions for the lower 
bounds of v* for the interference adaptive DC A schemes. 
However, computation of ;/* lin through these propositions 
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Fig. 8. Comparison of deterministic downlink lower bound with 
Zander-Ericksson' s result. 


is nontrivial due to the calculation of a{hk) j a s ihk) ■ A 
closed-form solution for a(hk) and a s (hk) is still a topic of 
future study. The results shown in this section are obtained 
numerically. All of bounds are computed for 7 = 18 dB and 
a = 4.0. In the shadowing model, a is assumed to be 8 dB. 

To compare with the lower bound of reuse partitioning DCA 
and v* of FCA proposed by Zander and Eriksson [20] (with 
channel reuse factor N — 9). 5 Fig. 8 depicts the downlink 
lower bound for v* . Two approximate bounds obtained from 
Theorem 2 are also shown. The lower bound is calculated from 
Theorem 3. The aihk) is obtained by interpolation. It is fairly 
understandable that the two bounds are very close when p is low, 
since in [20] the assumption of continuous reuse partitioning is 
used, which is impractical in location adaptive DCA, but is im- 
plemented by interference adaptive DCA. For a higher value of 
p, the Zander-Eriksson's bound is too loose in the sense that v* 
goes up to about 0.65 rather than 1 (as in proposed bound) when 
p goes to infinity. It should be pointed out that £ is pushed to in- 
finity in the definition of v* . As a result, the traffic dynamics 
among cells are eliminated. On the other hand, it is well known 
that FCA outperforms all traffic adaptive DCA schemes when 
£ is very large. Therefore, it is appropriate to view the differ- 
ence between the bound of IA-DCA and FCA in Fig. 8 as the 
gain from a more accurate representation of the CIR constraint 
in IA-DCA than in TA-DCA. 

The lower bounds of v* for the deterministic case and the 
shadowing case are displayed in Figs. 9 and 10, respectively. 
Comparison of the lower bounds in these figures shows that 
lognormal shadowing has some positive impact on the system 
capacity. For example, if we observe two downlink bounds at 
v* — 1%, the average traffic load p (user/ channel/ cell) of the 
system is approximately 0.36 in Fig. 10 and 0.29 in Fig. 9. This 
implies that 20% ~ 25% capacity gain can be obtained by uti- 
lization of shadowing with a typical value of variance (a — 8 
dB) for IA-DCA strategies. Even with two-way bounds, there 
is still about 8% ~ 10% capacity gain. This theoretical con- 
clusion is in agreement with the simulation results reported in 

5 7 =18 dB and a =4 .0 corresponds to hk « 7.22. Only the downlink 
bound in the deterministic case is available for comparison. 
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Fig. 9. Lower bound of v * for deterministic model. 
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n- 


hk h k £A k /A c 

hk (J2n-i A n )/A c = 1 implies that in order 

to handle the group of users with average number £, one needs 
up to K CRZ's with the largest reuse factor being Hk, From 
Fig. 11, it is known that the area of zones in ECRZ is always 
larger than the one in CRZ at the same hk ■ The value of hx de- 
creases from about 7.22 in the deterministic model to 5.56 in the 
shadowing case. This increase in Ak (or decrease in hx) results 
in the improvement of the bound for v* . 

Other conclusions that can be drawn from these figures 
include: 1) in the deterministic case, three bounds are very 
close to each other. This implies that interference on downlink 
and uplink channels are similar if channels are allocated 
by IAMP method and 2) in the shadowing case, the uplink 
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Fig. 11. Area change in CRZ and ECRZ. 

cochannel interference is the main limit to the capacity of the 
TDMA/FDMA system using IA-DCA schemes. 

VII. Conclusions 

In this paper, the computation of cochannel interference and 
outage probabilities is presented for both downlink and uplink 
channels in the TDMA/FDMA system using interference adap- 
tive DCA schemes. This computation is critical in the analysis of 
IA-DCA schemes. As an important application, the lower bound 
of the asymptotic probability of assignment failure is calculated 
for the downlink/uplink-only as well as two-way balanced in- 
terference adaptive DCA strategies. Both the deterministic and 
the shadowing propagation models are considered. The bound 
is shown to be a closed-form function of the traffic load p. The 
effects of CIR threshold 7 and the path-loss exponent a can 
only be observed in the numerical way. It is analytically shown 
that a capacity enhancement is possible by utilizing the shadow 
fading in DCA schemes. This capacity improvement cannot be 
achieved by either traffic adaptive or location adaptive DCA 
schemes, while it can be realized by interference adaptive DCA 
schemes. 

Our research plan is to find the performance bound of a real 
TDMA/FDMA system with the IA-DCA strategy by alleviating 
the assumption of arbitrarily large M, Ongoing work also in- 
cludes the application of the results in this paper to be used for 
more realistic traffic models, such as nonuniform traffic distri- 
bution. 
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